Isospin Diffusion in Heavy-Ion Collisions and the Neutron Skin Thickness of Lead 
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The correlation between the thickness of the neutron skin in 208 Pb, and the degree of isospin 
diffusion in heavy-ion collisions is examined. The same equation of state is used to compute the 
degree of isospin diffusion in an isospin-depedent transport model and the neutron skin thickness in 
the Hartree-Fock approximation. We find that skin thicknesses less than 0.15 fm are excluded by 
the isospin diffusion data. 
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I. INTRODUCTION 

The nuclear symmetry energy and its dependence on 
density holds a unique place in nuclear physics and as- 
trophysics. It is an essential ingredient in understanding 
many aspects of nuclear physics, from heavy- ion collisions 
to nuclear structure, and astrophysics, from supernovae 
and neutron stars to r-process nucleosynthesis (for a re- 
cent review see Ref. pj). At the same time, the nuclear 
symmetry energy is relatively uncertain; it is known to 
within only about 5 MeV at saturation density (com- 
pared to the binding energy which is known to within 
1 MeV). The density dependence of the nuclear symme- 
try energy at both sub-saturation and super-saturation 
densities is very poorly known. 

The symmetry energy is also an important quantity 
in determining the compressibility of matter. The in- 
comprcssibility K of symmetric nuclear matter at its 
saturation density p = 0.16 fm -3 has been determined 
to be 231 ± 5 MeV from nuclear giant monopole reso- 
nances and the equation of state (EOS) at densities 
of 2po < p < 5po has been constrained by measurements 
of collective flows 0, 0, B S m relativistic heavy-ion 
collisions §. However, further progress in determining 
more precisely both the parameter K and the EOS of 
symmetric nuclear matter is severely hindered by the un- 
certainties of the symmetry energy 0, . 

The possibility of determining the EOS of nuclear mat- 
ter from heavy-ion collisions has been discussed for al- 
most 30 years, and indeed impressive progress has been 
achieved. A number of heavy-ion collision probes of the 
symmetry energy have been proposed includin g is ospin 
fractionation [3 IH IH EE 111, isoscaling [lf| U^ . 
neutron-proton differential collective flow |l7j . pion pro- 
duction [l8j and neutron- proton correlation function [Ig ■ 
Determination of the equation of state from heavy-ion 
data is typically performed thr oug h comparisons with 
transport model simulations 0, |2(J, 0, [23, H3|. A re- 
cent analysis of the symmetry energy comes from isospin 
diffusion data for reactions involving 112 Sn and 124 Sn 
from the NSCL/MSU 0, |2^. Using an isospin- and 
momentum-dependent transport model, IBUU04 [26H27j . 
the NSCL/MSU data on isospin diffusion were found 



to be consistent with a nuclear symmetry energy of 
E S ym(p) ~ 31.6(p//5o) 105 at sub-saturation densities. 



The pressure of neutron-rich matter at densities near 
two-thirds saturation density (which is closely related to 
the symmetry energy, see Ref. Q) is tightly correlated 
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to the neutron skin thickness, SR, in 208 Pb [28ll29| . The 
present theoretical uncertainty in the symmetry energy 
results in a large uncertainty in the theoretical calcula- 
tions of the skin thickness. Experimental measurements 
of the skin thickness are also very uncertain, mostly re- 
sulting from the difficulty in disentangling the strong in- 
teractions between the probe and the neutrons in the 
nucleus. Values from 0.07 fm to 0.24 fm have been sug- 
gested by the data jU OOI • A determination of the sym- 
metry energy at sub-saturation densities would offer a 
prediction of skin thickness of lead or vice versa. 



The structure of nuclei and neutron stars and the prop- 
erties of heavy-ion collisions are all determined by the 
same underlying EOS. In particular, both the size of the 
neutron skin and the degree of isospin diffusion in heavy- 
ion collisions at intermediate energies are sensitive to the 
symmetry energy at sub-saturation densities. In this ar- 
ticle, we study the correlation between the thickness of 
neutron skin in 208 Pb and the strength of isospin diffusion 
using the same set of equations of state for asymmetric 
nuclear matter. This differs from previous work in that, 
for the first time, we make a direct connection between 
observables from nuclear structure and heavy-ion colli- 
sions. The transport model IBUU04 was used to predict 
the isospin diffusion for several equations of state which 
differ only in the density dependence of the symmetry en- 
ergy. These equations of state are also used to calculate 
5R for 208 Pb in the Hartree-Fock (HF) approximation. 
The correlation between the SR and the isospin diffusion 
provides a more stringent constraint than their individual 
values on the EOS of neutron-rich matter. Comparisons 
with the available data are also discussed. 
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II. THE EQUATION OF STATE AND THE 
TRANSPORT MODEL CALCULATIONS OF 
ISOSPIN DIFFUSION IN HEAVY-ION 
COLLISIONS 

Isospin diffusion measures quantitatively the net ex- 
change of isospin contents between the projectile and tar- 
get nuclei. Using symmetric reactions A + A and B + B 
as references the degree of isospin diffusion in the asym- 
metric reaction of A + B can be measured by [32^ 

_ 2X A+B - X A+A - X B+B 

Ri - x A+A - X B+B ^ ' 

where X is any isospin-sensitive tracer. In the recent 
NSCL/MSU experiments with 124 Sn on 112 Sn at a beam 
energy of 50 MeV/nucleon and an impact parameter 
about 6 fm, the isospin asymmetry of t he p rojectile-like 
residue was used as the isospin tracer j24|. Consistent 
with the experimental selection, in model analyses the av- 
erage isospin asymmetry (5) of the Sn-like residue was 
calculated from nuclcons with local densities higher than 
po/20 and velocities larger than 1/2 the beam velocity in 
the cm. frame. Reactions at intermediate energies are 
always complicated by preequilibrium particle emission 
and the production of possibly neutron-rich fragments at 
mid-rapidity, however, the quantity Ri has the advantage 
of minimizing significantly these effects [24| . 

The NSCL/MSU data were recently analyzed by using 
the IBUU04 version of an isospin and momentum depen- 
dent transport model [26ll2^|. Here, we recall the ingre- 
dients of the model that are most relevant for the present 
study. In this model one can select to use either the ex- 
perimental nucleon-nucleon cross sections in free space or 
reduced in-medium nucleon-nucleon cross sections. Since 
the later involves several model-dependent assumptions, 
we have used the free-space cross sections (as in Ref. pjjh 
The EOS used in this model is based on the Gogny effec- 
tive interactions. The potential U (p, S, p, r) for a nucleon 
with isospin t (1/2 for neutrons and —1/2 for protons) 
and momentum p in asymmetric nuclear matter at total 
density p is given by Hi HHH HI 

Umdi(p,S, P ,t) = A u ^ + aA + b(^-) (1-xS 2 ) 
Po Po \PoJ 
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where p T and p' T denote proton or neutron density with 
t / r'; and 5 = (p n — p p )/p is the isospin asymmetry. 
The / r (r, p) denotes the phase-space distribution func- 
tion at coordinate r and momentum p. The correspond- 
ing momentum-dependent interaction (MDI) leads to an 
incomprcssibility of Kq = 211 McV for the symmetric 



nuclear matter at saturation density. On the right hand 
side of Eq.©, the first four terms with cr = 3/4 and 
B = 106.35 MeV describe the momentum-independent 
interaction. The terms with parameters C r , r = —11.7 
MeV and C r , T ' = —103.4 MeV describe the momentum- 
dependent interaction of a nucleon of isospin r and mo- 
mentum p with like and unlike nucleons in the back- 
ground fields, respectively. We stress that the mo- 
mentum dependence in both the isoscalar and isovec- 
tor potentials is taken into account in a way consistent 
with that known empirically from nucleon optical po- 
tentials jl^l ■ With the parameter A = 1.0p F , where p F 
denotes nucleon Fermi momentum at po, the isoscalar po- 
tential (U n (p,p) + U p (p,p))/2 coincides with predictions 
from the variational many-body theory usin g inp uts con- 
strained by nucleon-nucleon scattering data |37ll38l|. and 
the isovector potential (U n (p,p) — U p (p,p))/2 also agrees 
with the momentum dependence of the Lane potential ex- 
tracted from nucleon-nucleus scattering experiments up 
to about 100 MeV [3!| an d (p ^n) charge exchange reac- 
tions up to about 45 MeV [37ll40l| . 

The parameter x in Eq. @ is introduced to allow vari- 
ations in the density dependence of the nuclear symme- 
try energy E syul (p), which is defined via the parabolic 
approximation to the nucleon specific energy in isospin 
asymmetric nuclear matter |2fil2lj. i.e., 

E(p, S) = E(p, 6 = 0) + E sym (p)6 2 + 0(c5 4 ). (3) 

The symmetry energies for three of the models consid- 
ered here are displayed in Fig.^ a s a function of density. 
With x = 1, E sym (p) matches what is predicted by a 
Hartree-Fock calculation using the Gogny effective inter- 
action H3 ■ The parameters Ai (x) and A u (x) arc Ai (x) = 
-120.57+2 J Ba;/(cr+l)and A u {x) = -95.98-2Bx/{a+l), 
respectively. It is important to note that large positive 
(negative) values of x lead to a "higher" ( "lower" ) sym- 
metry energy at sub-saturation densities, but a lower 
(higher) symmetry energy at super-saturation densities. 

An example of the results from the transport calcula- 
tions, taken from Ref. j2|| is given in Fig. for the MDI 
EOS with x = — 1. The time evolution of Ri is given for 
this equation of state. The evolution of the average cen- 
tral density p in units of the saturation density po is also 
given. It is clear that the late-time evolution of Ri for the 
MDI equation of state matches the data obtained from 
MSU. The central density at these times is quite small, 
indicating clearly that Ri is sensitive to the symmetry 
energy at sub-saturation densities. 

III. THE EOS FITTING PROCEDURE AND 
PROPERTIES OF 208 PB 

Because the MDI EOS in the form above does not eas- 
ily avail itself to the calculation of the structure of lead 
nuclei we choose to fit the EOS to a Skyrme model and 
use a Skyrme-Hartree-Fock code to calculate correspond- 
ing lead nucleus. This procedure has been used success- 
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FIG. 1: The symmetry energies as a function of baryon num- 
ber density, p, in units of the saturation density, po, for three 
of the equations of state considered in this work. 
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FIG. 2: Time evolutions of Rt and average central density for 
MDI and SBKD interactions with x = — 1. 
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FIG. 3: The results of fitting the MDI EOS for x = -1 (solid 
lines) to a Skyrme EOS (dashed lines). Plotted is the energy 
per baryon as a function of the baryon density, p. The dashed 
lines are barely visible because of the high quality of the fit. 

x E b (MeV) R ch (fm) SR (fm) L (MeV) K aym (MeV) 
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62.1 
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0.16 
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Exp. 
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0.07-0.24 







TABLE I: The properties of lead nuclei for the Skyrme equa- 
tions of state and their experimental values. L and K 3ym 
are the slope and curvature of the corresponding symmetry 
energy at po. 



fully to predict the neutron skin based on the Akmal, 
et. al. EOS H3 in Ref. Q. The four MDI equations of 
state with different values of x were fit to Skyrme mod- 
els constrained so values for the binding energy (Eb) and 
charge radius (R c h) for lead which were within 2% of 
the experimental values. For the equation of state with 
x = — 1, the results of the fit are displayed in Fig. 
The similar nature of the Skyrme and MDI parameteri- 
zations allows for a very close fit. For all of the four fits 
employed, the slope, L = 3po (dE sym / dp) , and the 
curvature K sym = 9pp {d 2 ''E synl / dp 2 ) of the symme- 
try energy and the final properties of lead nuclei in the 
Hartree-Fock approximation are given in Table [Q Also 
given in the table are the experimental values. 

The effective mass for the original MDI EOS and the 
new Skyrme model are not as close, for example, the ef- 
fectve mass at saturation density for the MDI EOS with 
x = — 1 is 0.67, while for the Skyrme fit, it is 0.77. How- 
ever, the neutron skin thickness is more sensitive to the 
bulk energetics than the effective mass. As a demonstra- 



tion, a new Skyrme fit to the MDI EOS with x = — 1 
with a lower effective mass of 0.65 only lowers the skin 
thickness by 0.01 fm. 



IV. ISOSPIN DIFFUSION VERSUS THE 
NEUTRON-SKIN OF 208 PB 

Now we turn to the correlation between the degree of 
isospin diffusion and the size of neutron-skin in 208 Pb. 
Do we expect a correlation between these two seemingly 
different observables? As it has been discussed in detail 
in Refs. |42ll43l| . the size of neutron skins in heavy nuclei 
is determined by the difference in pressures on neutrons 
and protons 8P. The latter is propotional to L, which 
measures the stiffness of the symmetry enrgy at satura- 
tion density. As shown in Table [IJ the symmetry energy 
becomes more stiff (for densities larger than the satura- 
tion density) as x decreases from 1 to —2, and the size 
of the neutron skin in 208 Pb increases. Alternatively, one 
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can understand this effect by noting that the the symme- 
try energy is lower (e.g., with x = —2) at subsaturation 
densities the cost of creating a difference between the 
neutron and proton densities is small and the skin thick- 
ness is large. As the symmetry energy increases with 
larger values of x, the neutron skin thickness decreases 
accordingly. 

It is also well known that the degree of isospin diffusion 
in heavy- ion collisions depends sensitively on the stiffness 
of the symmetry energy p2l I24I I2MI44I KBj . A correlation 
between the degree of isospin diffusion in heavy-ion colli- 
sions and the size of neutron skins in heavy nuclei using 
the same EOS is thus expected. Of course, it is under- 
stood that this correlation is not universal as both the 
isospin diffusion and size of neutron skins are system de- 
pendent. However, an examination of this correlation for 
any system and simultaneous comparisons with the cor- 
responding data may provide a more stringent constraint 
on the underlying EOS. Shown in Fig.0|is our analysis of 
this correlation for the isospin diffusion in 124 Sn+ 112 Sn 
and the neutron skin in 208 Pb. The available data on 
both quantities are also included in the graph. There are 
some remaining systematic uncertainties that are difficult 
to estimate which are not addressed here (see Ref. j2(| 
for more details). 

From Eq. we can see that R4 = 1 implies the 
projectile-like residue has the same isospin asymmetry as 
the projectile 124 Sn without any net exchange of isospin 
asymmetry with the target 112 Sn. This is also often re- 
ferred as the complete isospin transparency. A value of 
zero for Ri implies that a complete mixing between the 
projectile and target indicating the establishment of an 
isospin equilibrium during the reaction. 

The NSCL/MSU data indicates that the mixing is 
about 50%. Unlike the relation between x and the neu- 
tron skin thickness, the relation between Ri and the pa- 
rameter x is non-monotonic. This is because the Ri is di- 
rectly related to the symmetry potential which depends 
on both the density and momentum as discussed in detail 
in Ref. [25J. Transport model calculations with x = — 1 
are closest to the experimental data, although x = and 
x = — 2 are not necessarily excluded within the current 
uncertainties of both the experimental data and model 
calculations. 

Upon comparing with previous estimates for 5R, the 
connection with the moderate amount of isospin diffusion 
displayed in heavy-ion reactions suggests that SR is on 
the upper range of the present experimental uncertainty. 
Only skin thicknesses larger than about 0.2 fm are con- 
sistent with the transport model calculations. Skin thick- 
nesses lower than 0.15 fm appear to be excluded. How- 
ever, the isospin dependence of the in-medium nucleon- 
nucleon cross sections may affect the isospin diffusion. It 
was shown analytically that the degree of isospin diffusion 
depends on both the symmetry e nerg y and the neutron- 
proton scattering cross section j44|- Because of the 
momentum-dependence of the single nuclcon potential of 
eq.IH not only the values of the n-n and p-p cross sections 
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FIG. 4: The skin thickness in lead, 5R, versus the isospin 
diffusion parameter 1 — R4 for the four equations of state 
in this work. A sample measurement with an uncertainty 
corresponding to a 1% measurement of the neutron radius 
is included, as well as the present acceptable experimental 
range. 



are reduced and different from each other compared to 
their free-space values, the ratio 2o~ np / (o~ pp + a nn ) is also 
modified from their free-space value. Whic the reduc- 
tion of in-medium n-n and p-p scattering cross sections 
is expected to have little effect on the value of Ri, the 
reduced in-medium n-p scattering cross sections leads to 
smaller isospin diffusion. To reproduce the NSCL/MSU 
data a larger symmetry energy at sub-saturation densi- 
ties is then needed. The use of in-medium cross sections 
will, in effect, decrease the preferred neutron skin thick- 
nesses for the isospin diffusion data. An investigation of 
this effect is currently underway and the results will be 
reported elsewhere. 

Although isospin diffusion in heavy-ion collisions pri- 
marily constrains the symmetry energy at sub-saturation 
densities, it is useful to check to ensure that the infor- 
mation from heavy-ion collisions is not inconsistent with 
neutron star observations. Neutron star radii are sensi- 
tive probes of the symmetry energy at higher density |46j . 
We solve the Tolman-Oppcnhcimcr-Volkov equations for 
the MDI EOS with x = 0, -1, and -2. The maximum 
masses (in order of decreasing x) are 1.92, 2.08, and 2.13 
Mq and the radii of the maximum mass stars are 10.1, 
11.2, and 11.8 km. In addition to being consistent with 
the compactness constraint given by the measurement of 
the redshift of EXO0748-676 in Ref. H3, this equation 
of state is also consistent with the recent neutron star 
radius measurement in Ref. |48[ . 

Recently, Ref. nas a l so explored the impact of the 
symmetry energy on the neutron skin thickness in lead. 
Using an relativistic mean-field model they construct 
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an equation of state which faithfully describes the giant 
monopole and dipole resonances in nuclei. Their model 
suggests a neutron skin thickness of 0.21 fm in Pb 208 
which is consistent with our constraints from isospin dif- 
fusion in heavy-ion collisions. 

V. SUMMARY 

In summary, we examined the correlation of the de- 
gree of isospin diffusion in heavy-ion collisions and the 
size of neutron skin in 208 Pb using the same EOS within 



the IBUU04 transport model and the HF approach. We 
found that neutron skin thicknesses less than 0.15 are 
disfavored by the isospin diffusion data. Reduction of 
the experimental and theoretical uncertainties will likely 
lead to an important constraint on the symmetry energy 
at sub-saturation density. 
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